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Abstract 

The surface stabihty of two interacting (for example, by van der Waals forces) incom- 
pressible thin films, one bonded to a substrate and the other to a contactor, is studied 
extending the work of Shenoy and Sharma, Physical Review Letters 18, 119-122 (2001). 
The analysis indicates that the wavelength of the instability depends strongly on the shear 
moduli and thicknesses of the films but not on the nature and magnitude of the interaction. 
When the films have equal shear moduli, the wavelength of the instability has an intermedi- 
ate value between the wavelengths of the instabilities had each of the films been interacting 
with rigid contactors. On the other hand, if the films have different shear moduli but equal 
thicknesses, then the wavelength of the instability is identical to that had the films been 
interacting with rigid contactors. In the more general case when the two films have different 
shear moduli and thicknesses, the nature of the critical wavelength is more complex. When 
ratio of the shear moduli of the contacting film to that of the film bonded to the substrate 
exceeds 5.32, the wavelength of the of the instability jumps from the value close to that 
determined by the thickness of the film bonded to the substrate to that of the contacting 
film, as the thickness of the contacting film is increased. 

1 Introduction 

The mechanics of adhesion and contact between two elastic bodies has attracted attention of 
researchers over the years. Following the classical work of Hertz, an important step in this 
area was taken by Johnson, Kendall and Roberts (1971) who identified the importance of the 
interactions between the contacting surfaces. These interactions can be due to van der Waals 
forces, electrostatic forces etc. between the contacting surfaces. Indeed, measurements of van 
der Waals forces between surfaces (Tabor and Winterton 1969) were available when Johnson et 
al. (1971) formulated their famous JKR theory of contact. 
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More recently, two groups (Ghatak, Chaudhury, Shenoy and Sharma 2000, Monch and Her- 
minghaus 2001) have reported experiments aimed at understanding contact and adhesion me- 
chanics between two elastic bodies with planar topology. The experiments reported by Ghatak et 
al. (2000) used a configuration where the glass plate was placed in contact with the elastomeric 
film between two spacer bars creating a small gap between the film and the glass plate. Monch 
and Herminghaus (2001) performed experiments that involved the contact of a glass plate with 
an elastomeric thin film bonded rigidly to a glass substrate. In both these experimental works it 
was observed that the film lost planarity when the contacting glass plate reached contact prox- 
imity (10-50nm), resulting in a pattern with a well defined wavelength. The two key features 
observed in both sets of experiments are (a) the wavelength of the pattern depends linearly on 
the thickness of the film (b) the wavelength is not affected by the magnitude and nature of the 
interactions. 

The theoretical interpretation of this results were first reported by Shenoy and Sharma (20016) 
who argued that the instability occurs due to a competition between the interaction energy of 
the film with the contactor and the elastic energy in the film due to inhomogeneous deformation. 
They showed that instability sets in when the ratio of the "stiffness of interaction" and the elastic 
stiffness (defined as shear modulus divided by the thickness) exceeds a critical value. In addition, 
they pointed out that the wavelength of the instability is determined solely by the elastic energy 
of the film and explained the linear dependence of the wavelength on the thickness of the film. 
A detailed account of this work can be found in Shenoy and Sharma (2001a). The origins and 
nature of this instability are different from those known to exist in solid films (Asaro and Tiller 
1972, Grinfeld 1993, Srolovitz 1989, Ramirez 1989) and fiuid films (Herminghaus, Jacobs, Mecke, 
Bischof, Fery, Ibn-Elhaj and Schlagowski 1998, Sharma and Khanna 1998, Reiter, Khanna and 
Sharma 2000, Schaffer, Thurn-Albrecht, Russel and Steiner 2000). 

The purpose of this paper is to extend the work of Shenoy and Sharma (20016) to the case 
when the contactor also has a film bonded to it as shown in fig. 1. The contacting film may 
have different properties in that the shear modulus and thickness may be different from the film 
bonded to the substrate. It is shown that one of key features i. e., the independence of the 
wavelength of the instability on the nature of interaction is unchanged. The critical wavelength, 
although does not depend on the nature and strength of the interaction, does depend strongly 
on the thicknesses and moduli of the film. It is shown that the relative stiffnesses of the films 
play an important role in the determination of the wavelength of the instability when the films 
have widely differing properties; the complete dependence of the wavelength on the thicknesses 
and moduli of the film is obtained. This paper is written with a purpose of motivating further 
experiments in the kinds of systems indicated in fig. 1 and is organised as follows. The next 
section contains a stability analysis of two interacting fihns. Results of the analysis are presented 
and discussed in section 3. The findings of this paper are summarised in section 4. The appendix 
at the end of the paper contains a short calculation of the normal traction along the surface of 
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an incompressible film subjected to a sinusoidal deformation. 

2 Stability of Two Interacting Films 
2.1 Model Description 

A schematic depiction of model considered for the study of stability of interacting thin films is 
shown fig. 1. The system consists of a substrate-contactor configuration where both the substrate 
and the contactor have thin films (possibly of different materials) bonded to them. The substrate 
film (called film a) has a thickness ha, while the contactor film (called film b) has a thickness hi,. 
A coordinate system described by coordinates {xi,X2) is used to describe position vectors. The 
positive '2' direction is the outward normal of film a, while the negative '2' direction point in 
the normal direction to the surface of the film b. We consider only plane strain deformation of 
the system. 

The equilibrium configuration of this system is determined by the potential energy 



where W is the strain energy density, e is the strain tensor, u is the displacement vector (the 
superscript a denotes value in film a etc.) with V being an appropriate measure of the volume 
(includes volumes of both the films), and S is the intcrfacial area of the two films over which they 
interact. The function U represents the interaction potential between the two films; it is this 
term that gives rise to interesting physics in this system. We make two key physical assumptions 
in writing the total potential energy: (i) the contribution from the surface energies of the films 
is negligible (ii) the films are considered to be made of incompressible elastic materials. Both of 
these assumptions are valid in physical systems where such instabilities are triggered, as shown 
by (Shenoy and Sharma 20016). 

The physical origins of such an interaction can be any of the several causes - van der Waals 
forces, electrostatic forces between films etc; the potential t/ is a generic interaction potential. If 
the potential is due to the attractive van der Waals interaction, then U is described by 



where A is the Hamaker constant (of the order of 10~^^J), and do is the distance between the 
surfaces of the two films, i. e., the gap thickness. The strength of the interaction potential is 
determined by the gap thickness dg- The contactor is imagined to be brought in the proximity 
of the substrate by reducing the distance do - the interesting quantity being the distance dc at 
which the attractive interactions are strong enough to trigger instability in the system. 

Linear stability analysis is performed using linear kinematics and a linearised interaction 
potential. To this end, the interaction potential is expanded in a power series about the reference 
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Figure 1: A thin elastic film bonded to a rigid substrate interacting with another film bonded on 
to a contactor. The dashed lines shows possible inhomogeneous deformation of the films when 
instability sets in. (Distances are not to scale.) 

state of the undeformed films and terms of up to quadratic order in (w" — u^) ■ n are retained 



where 



U{{u'' -u^)-n)^Uo + F (u" -u^)-n + -Y ((w" - u^) ■ nf 
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t/o = f/(0), F = f/'(0), Y = U"{Q). 



(3) 



(4) 



The quantity F, called the interaction stiffness, is of importance and governs the stability of the 
system. The aim of the analysis is to find the condition(s) on Y under which instability sets in. 
The above approximations give an expression for the approximate potential energy as 



n„ = y W{e)dV -J (^Uo + F {u" -u'')-n + ^Y (K - u^) ■ 
The strain energy density W is 

in film a 



(5) 
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in film b 



where fia and fib are respectively the shear moduli of the films a and b. The equilibrium displace- 
ment fields in the films minimise the potential energy (5) while satisfying the rigid boundary 
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conditions at the film-substrate interface in film a, and the film contactor interface in film b. In 
addition, the stresses derived from these displacement fields satisfy the condition of vanishing 
shear stress ai2 at the surface in both the films. The normal stresses satisfy the condition 

<^22 = F + Y{u^,-4) (7) 



4, = F + Y{u^,-ul) (8) 

at their respective surfaces. 

Homogeneous Solution: Since both the films are incompressible, the homogeneous solution has 
displacements vanishing everywhere in both the films, and stress state in both the films is one 
of constant pressure (equal to F). The point in question is the stability of this homogeneous 
solution, i. e., for what value(s) of Y does instability occur? 

2.2 Stability Analysis 

To study the stability of the system above, the homogeneous solution is perturbed by bifurcation 
fields (denoted by Ui here and henceforth) such that the surfaces of the film a and b, respectively, 
have displacements 

■^2(^1) — acos{kxi) (9) 



M^(a;i) = Pcos{kxi). (10) 
The additional stresses produced by these fields satisfy 

<2 = 0, = 0, (11) 

and 

^22 = y{y'^2 - 4) (12) 



<^22 = ^« - 4) (13) 

along the interacting surfaces of the films. All symbols have obvious meanings. 

It is shown in the appendix that the stress (T22 the surface of a film a whose displacement 
is U2{xi) = a cos(A;a;i) with 17^2 = is given by 

^22(^1) ^'2^"' S{hak)kacos{kxi), (14) 

and, the stress (722(3^1) along the surface of the film b is 

a^^ixi) = -2p!'S{hbk) kf3cos{kxi), (15) 
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where 5" is a nondimensional function defined as 

l + cosh(2e) + 2e' 
= smh(20-2e ■ f^"* 

For future discussion we note that S'(^) — > 1/^^ as ^ — > 0, and S{^) — > coth(^) when ^ — > oo. 

Substitution of (14) and (15) in (12) and (13) leads to the following homogeneous set of 
equations for a and (3: 

'2pLakS{hak)-Y Y = r° 

-Y Y -2f,t,kS{hhk) J \P J [O 

The condition for the existence of nontrivial bifurcation fields, i. e., the condition for the existence 
of nontrivial solutions for a and (3 is obtained by setting the determinant of the coefficient matrix 
to zero resulting in a relation between the interaction stiffness Y and the wavenumber k of the 
bifurcation field 

2kfXafXbS{hak) S{hbk) 

HaS{hak) + ljLbS{hbk) ' 

If, for a given value of Y, there is at least one real value of k that solves (18), then the ho- 
mogeneous solution is unstable and the films deform inhomogeneously. The lowest value of Y 
for which bifurcations are possible is called the critical interaction stiffness denoted by Y^, the 
wavenumber(s) of the mode(s) that satisfies (satisfy) (18) for Y — Yc is(are) called the critical 
mode(s) and the(these) wavenumber(s) is(are) denoted by kc- 

At this point, there are four parameters /la, A*6, ha, hb that enter into the determination of the 
stabihty of the system. The discussion to follow is much simphfied on introduction of effective 
parameters; to this end, we define 

/x= — - — 19 

fJ'a + fJ'b 

and 

h^K + hb. (20) 
We also introduce nondimensional parameters M and H defined as 

A = M, = (1 - M) (21) 

and 

ha = Hh, hb^{l- H)h. (22) 

When Ha <^ pib: M ^ 1 and M ^ if ^ Ha', similarly when ha ^ hb, H ^ 1, while if — >■ 
implies hb ha- Another important quantity of interest is the effective elastic stiffness K^ff of 
the two-film system defined as 

hah ^ ^a^b ^fi 1 , 

-ff /^^A^ i^ahb + f^bha h{{l-H){l-M) + HMy ^ ' 

ha hb 
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Based on the above definitions, the expression for the interaction stiffness in (18) can be 
recast in nondimensional form as 

Y _ 2 ((1 - - M) + HM)qS{Hq) S{{1 - H)q) 



(24) 



immediately recovered. When M — > 1, — > ^ and 

ha 



Keff {l-M)S{Hq) + MS{{l-H)q) 

where q — hk. 

The results for the case when //^ » obtained by Shenoy and Sharma (2001 &) can be 

n M 

^^2Hq S{Hq) = 2 hak S{hak) (25) 

The critical interaction stiffness can be obtained as Yc = 6.22fia/ha and hake = 2.12, which is 

precisely the result of Shenoy and Sharma (20016). Another route to the same result is to let 

ha ^ hb, i. e., as H ^ 1. In this case, i^e — ^ t~- again, and 

ha 

^ = 2Hq S{Hq) = 2 hak S{hak) (26) 

where the property of the function S near 0, noted above, is used. It can be easily shown that 
similar results are obtained when either //^ <^ Ha {M 0) or hi, ^ ha (i? — > 0). 

In this section, we have derived the condition for the onset of instability in a system with 
two interacting incompressible elastic films. In addition, it is shown that the results of Shenoy 
and Sharma (20016) who treated a film interacting with a rigid contactor, are a special case of 
the present formulation. The more general case is taken up for discussion in the next section. 

3 Results and Discussion 

The nondimensional quantities M and H introduced in the last section can take any value from 
to 1. However, it is evident that the results for a given value of M and H are physically identical 
to that of (1 — M) and (1 — H) (as is evident from (24)). Due to this reason, only the regime 
I < M < 1 and < if < 1 is considered. Results are discussed in three categories: (i) //„ = A*6, 
ha 7^ hb (ii) /la 7^ A*6, ha — hb (iii) i^a 7^ A*6, ha 7^ hb- The results of all these cases are plotted in 
two sets of graphs shown in fig. 2 and fig. 3. 

An important point about the results is that, just as in the case of a film interacting with a 
rigid contactor, the wavelength of the instability does not depend on the magnitude and nature 
of the interaction forces, but depends only on the thicknesses of the films and their moduli, i. e., 
on the parameters H and M. 

3.1 Films with Equal Shear Moduli 

When the shear moduh of the two films are equal (/^a = jjib), M = \^ ji — jia/'^ and K^ff — 
H Hal ha- The results for Yc and kc correspond to M = 0.5 in fig. 2 and fig. 3. The magnitude 
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Figure 2: The critical interaction stiffness as a function of H for various values of M. 
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of non dimensional critical interaction stiffness Yc/Kf.ff increases with decreasing H attains a 
peak and falls as H approaches |. If /la is kept fixed and hf, is changed to change if, then 
increasing Yc/K^ff does not imply an increase in Y^, since -f^e// f^-Us linearly with H. Thus 
for a fixed ha, the effective stiffness of the system reduces with decreasing H. For example, 
when H — 3/4, i.e, hb/ha — 1/3, the value of — 3ijLa/4:ha and the value of Yj — 6.92 or 
Yc — B.lQ/ia/ha which is less than Yc — 6.22//a//ia when the film a interacts with a rigid contactor. 
Looking at the same example from the point of view of film b, we note that Ke — {1 — H)iih/hb 
and Yc — l.lS/ib/hb which is again less than the value 6.22iih/hb. In fact, when films have 
equal thicknesses and equal shear moduh, the value of Yc — SAlHa/ha. The conclusion of these 
observations is that the magnitude of the critical interaction stiffness is less than that when either 
of the film is interacting with a rigid contactor. Physically, this can be understood by the fact 
that introduction of a second film on the contactor makes the system more compliant, i. e., the 
effective elastic stiffness of the system comes down. Thus, the elastic energy penalty required 
to cause an inhomogeneous deformation in the system also comes down. As a consequence, the 
absolute value of the interaction stiffness required to trigger the instability is reduced. Moreover, 
the actual value of Yc is governed by the more compliant film. 

Turning now to a discussion of the wavenumber of the critical mode, we note that as H 
decreases from 1, hake increases and then falls subsequently. For values of H close to 1, the 
instability triggers a shorter wavelength mode in film a as compared to when the film a interacts 
with a rigid contactor. For example, when H — 3/4, the value of hake = 2.45; in terms the film b 
the value of kc is hbkc = 0.82, i. e., the wavelength of the critical mode is much larger than that 
if the film b had been interacting with a rigid contactor. The key point is that the wavelength of 
the instability in this regime is intermediate to critical wavelengths determined by the thicknesses 
of the participating films. As the thicknesses are made equal, i. e., when i7 = |, the critical 
wavelength becomes hake — 2.12 which is exactly the wavelength had either of the films been 
interacting with a rigid contactor. When M = 1/2, the case of H < 1/2 is physically identical 
to the case H > 1/2. These results can be interpreted physically as follows. The dependence 
elastic energy in the film on the wavenumber k is determined solely be the thickness of the film, 
more precisely by the function S (defined in (16)) which depends on the thickness of the film, 
the shear modulus is a mere multiplicative parameter. The elastic energy (per unit length of the 
film) is a minimum when the wavenumber is equal to 2. 12/ (thickness). However, when there are 
two films with very different thicknesses, the total elastic energy in the system (sum of elastic 
energies in the two films) can be minimised by choosing an intermediate wavelength. 

3.2 Films with Equal Thicknesses 

When the films are made of equal thickness ha — hi,, the parameter H — \. For this case, the 
effective stiffness Ke — /i/ha, and the nondimensional critical interaction stiffness is independent 
of the value of M, with Ye — 6.22/1/ ha- Since /i < /la, and /i < /ib, the actual magnitude of 
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Figure 3: The wavelength of the critical mode as a function of H for various values of M. 
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the critical interaction stiffness is less than the critical interaction stiffness of the films if they 
had been interacting with rigid contactors. This can, again, be understood by the argument 
that the effective stiffness of the system comes down on introduction of the second film on the 
contactor. The critical wavenumber, however, does not depend on M as is evident from fig. 3 
and is equal to the wavelength had the films been interacting with a rigid contactor. Again, this 
can understood from the argument stated above that the elastic energy of a mode is determined 
by the thickness of the film via the function S and the moduli are multiplicative constants in the 
energy expression. This idea is evident from (24) on substitution oi H — ^. 

3.3 Films with Unequal Thicknesses and Moduli 

The physics of the instabihty in the general case can be understood based on the results of the 
previous two cases. When M > i, the film a is softer than the film 6, i. e., /Xq < A*&- When 
M is close to 1/2, the qualitative features of the instabihty are unchanged. The nondimensional 
interaction stiffness Yc/Keff increases with decreasing H from 1, and falls as H approaches |. 
On further decrease of H, Yc/K^ff increases and decreases again to 6.22 as H approaches 0. The 
behaviour of the critical wavelength is qualitatively identical to that of the case when M — \. 

As M is made close to 1 (M < 1), a quahtative change appears in the solution. In fact, 
when M — 0.8418 (corresponds to jibl l^-a — 5.32), the critical wavenumber changes suddenly as 
a function of H at about H — 0.1756 (corresponds to hi,/ha — 4.69). In fact, for larger values of 
M, as H is reduced, the critical wavenumber remains a constant (and close to the wavenumber 
hake — 2.12) but jumps suddenly to a much lower value - the value of H at which the jump 
occurs is called Hj. This is evident when M — 0.9 (shown in fig. 3) where the wavenumber of 
the critical mode hake ~ 1.79 when H < 0.25; however, at if = 0.095 = Hj (corresponding to 
hb/ha — 9.52), the critical wavelength becomes hake — 0.22 — > h^kc — 2.1, i.e., the wavenumber 
of the instability changes to the critical wavenumber of the film b had it been interacting with 
a rigid contactor. The value of Hj, i. e., the nondimensional value of H at which the critical 
wavenumber jumps from a higher value kh to a lower value ki depends on the ratio of the shear 
moduli; this dependence is plotted in fig. 5. It is evident that Hj decreases with M, as is expected 
from the physical arguments presented below. The two wavenumbers at Hj are plotted in fig. 6. 
Again, as M approaches 1, hakh approaches 2.12, the value determined by the thickness of the 
film a, and haki approaches (it can, in fact, be shown that hfjki approaches 2.12). 

Prom a physical perspective, the total elastic energy is the sum of elastic energies in the 
films which are products of their respective shear moduli with the S function defined in (16). 
When the elastic moduh are very different, the total energy of the system can have multiple 
minima, which implies that y as a function of q defined in (24) has multiple minima. Indeed, 
a plot of y as a function of hak shown in fig. 4 clearly shows the presence of multiple minima; 
additionally, the value of hak where the minimum is attained changes to a smaller value of hak as 
H is reduced. The key physical idea is that instabihty is governed by the film that has a smaller 
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Figure 4: The interaction stiffness Y/ K^ff as a function hak showing muhiple minima for different 
values of H. Note how the value of hak at which Y/K^ff attains a minimum changes as H is 
reduced. 
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Figure 6: The lower {hah) and higher (hakh) wavenumbers ect H — Hj, as a function of M. 

stiffness (the comphant film). Although /if, is much larger than hf, can be made larger than 
ha to the extent that the stiffness of the film b is smaller than that of a. Thus when H is lowered, 
the critical wavenumber of the instability of the system jumps from a wavenumber close to the 
critical wavenumber determined by film a (as if interacting with a rigid contactor) to that of the 
film b (as if interacting with a rigid substrate). 

4 Summary and Conclusion 

This paper extends the work of Shenoy and Sharma (20016), who studied the interaction of a 
soft thin film with a rigid contactor, to the case where the contactor also has a soft film bonded 
to it. This new analysis uncovers some interesting new results that are summarised below. 
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1. Instability sets in when the interaction stiffness exceeds a critical value which is determined 
by the shear moduh and thicknesses of the interacting films. The magnitude of the critical 
interaction stiffness is smaller than that of the critical interaction stiffnesses of the indi- 
vidual films had they been interacting with rigid contactors. Thus it is "easier" for the 
instability to occur in the two-film system. 

2. In all the cases, the wavelength of the instability does not depend of the nature and magni- 
tude of the interaction. The critical wavelength is determined only by the thicknesses and 
shear moduli of the films. 

3. When the films have equal shear moduli, the critical wavenumber of the instability depends 
on the relative thicknesses of the films. When the thicknesses of the films are very different, 
the wavelength of the instability is intermediate between the wavelength of the instability 
in each of the films had they been interacting with rigid contactors. However, when the 
thicknesses of the film are equal, the wavelength of the instability is exactly same as that 
when the films interact with rigid contactors. 

4. When the films have equal thicknesses, but different shear moduli, the magnitude of the 
critical interaction stiffness is less than that of each of the film had they been interacting 
with rigid contactors. The wavelength of the critical mode, however, is equal to that when 
the films interact with rigid contactors and is independent of the shear moduli of the films. 

5. When the films are of both unequal thicknesses and moduli, the behaviour is dominated 
by the less stiff films (stiffness being the ratio of the shear modulus to the thickness of 
the film) . As long as the modulus of the contacting film is less than about five times that 
of the substrate film, the behaviour of the solution is qualitatively same as that discussed 
in the above cases. However, when the modulus of the contacting film is larger than 5.32 
times the thickness of the film on the substrate, the wavelength of the instability depends 
strongly on the thicknesses of the films. As the thickness of the contacting film is increased, 
the wavelength of the instability jumps from a value determined by the thickness of the 
substrate film to that determined by the contacting film. 

It is hoped that the paper will stimulate further experiments along the lines of Ghatak et al. 
(2000) and Monch and Herminghaus (2001) to verify the results presented here. The analysis 
presented here neglects the effect of surface energies and the compressibility of the films. The 
inclusion of these effects along with the viscosity of the films will be presented elsewhere (Sarkar, 
Shenoy and Sharma 2001). 
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A Stresses Along the Surface of an Incompressible Film 
Bonded to a Substrate with Sinusoidal Surface Defor- 
mation 

The aim of this appendix is to outhne the determination of the normal traction along the surface 
of an incompressible film bonded to a rigid substrate, the surface of which is deformed sinusoidally. 
The main result of this appendix is the formula used in (14) and (15). 

Consider an incompressible elastic film described by coordinates (xi, ^2) of thickness h bonded 
to a rigid substrate such that the free surface of the film has coordinate ^2 = and the interface 
between the film and the substrate has coordinate X2 — —h. The film has a shear modulus /i. 
The boundary value problem has the following boundary conditions: 

ui{xi, -h) = U2{xi, -h) = 0, (27) 

U2{xi,0) — a cos{kxi), (7i2(xi, 0) = 0. (28) 

The equilibrium equation in terms of the displacements is 

II Uijj + p^i ^ 0, (29) 

where p is the pressure field. The incompressibility condition is expressed as 

Ui,i = 0. (30) 

A general solution of the set of differential equations that anticipates the boundary conditions 
(27) and (28) is 

ui(xi,X2) = ~{{B + k{A + Bx2))e'''=' + {D-k{C + Dx2))e-'''=')sm{kxi) (31) 
U2{xi,X2) = a{{A + Bx2)e''''^ + {C + Dx2)e-^''')cos{kxi) (32) 

p{xi,X2) = -2iia {Be''^^ + De-''''^) cos{kxi) (33) 

The constants A, B, C and D can be determined from the boundary conditions (27) and (28). 
The solution is 

1 + e^^^ - 2fcMl - fc/^) c-l-A (U) 
2smh{2kh)-Akh' /^-/c + i^. (35) 
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Prom these relations, the expression for the normal component of traction along the surface 
of the film can be derived as 



o'22{xi, 0) = 2/xM2,2 + P = 2ii{2A — l)ka cos{kxi). 



(36) 



On substitution of the solution for A from (34), the expression simphfies to 



(722 (2^1,0) = 2iiS{hk) ka cos{kxi) 



(37) 



where the function S is defined in (16). 
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